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METHOD, APPARATUS AND SYSTEM FOR REMOVING MOTION 
ARTIFACTS FROM MEASUREMENTS OF BODILY PARAMETERS 

PRIORITY CLAIM 
This application claims the benefit of the filing date of United States Provisional 
Patent Application Serial Number 60/104,422, filed October 15, 1998, for "METHOD 
FOR REMOVING MOTION ARTIFACTS FROM DEVICES FOR SENSING 
BODILY PARAMETERS AND APPARATUS FOR EFFECTING SAME". 

TECHNICAL FIELD 
This invention relates to the field of signal processing. More particularly, this 
invention relates to processing measured signals to remove unwanted signal 
components caused by noise and especially noise caused by motion artifacts. 



BACKGROUND 

The measurement of physiological signals can often be difficult because the 
underlying physiological processes may generate very low level signals. Furthermore, 
interfering noise is inherent in the body and the interface between the body and sensors 

20 of physiological processes. Examples of physiological measurements include: 

measurement of electrocardiogram (ECG) signals based on the electrical depolarization 
of the heart muscle, blood pressure, blood oxygen saturation, partial pressureof C0 2 , 
heart rate, respiration rate, and depth of anesthesia. ECG signals, for example, are 
typically detected by surface electrodes mounted on the chest of a patient ECG signals 

25 are weak at the signal source (i.e., the heart) and are even weaker at the surface of the 
chest. Furthermore, electrical interference from the activity of other muscles (eg., 
noise caused by patient breathing, general movement, etc.) causes additional 
interference with physiological signals such as an ECG. Thus, considerable care must 
be taken in the design and use of physiological processors to enhance the quality of the 

30 true signal and reduce the effects of interfering noise signals. 

It is convenient to characterize a measured signal as being a composite signal 
composed of a true signal component and a noise signal component. The terms 
"measured signal" and "composite signal" will be used interchangeably hereinafter. 
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Signal processors are frequently used to remove noise signal components from a 
composite measured signal in order to obtain a signal which closely, if not identically, 
represents the true signal. Conventional filtering techniques such as low pass, band 
pass, and high pass filtering can be used to remove noise signal components from the 
5 measured composite signal where the noise signal component occupies a frequency 
range outside the true signal component More sophisticated techniques for 
conventional noise filtering include multiple notch filters, which are suitable for use 
where the noise signal component exists at multiple, distinct frequencies, all outside the 
true signal frequency band 
10 However, it is often the case that the frequency spectrum of the true and noise 

signal components overlap and that the statistical properties of both signal components 
change with time. More importantly, there are many cases where little is known about 
the noise signal component In such cases, conventional filtering techniques are 
ineffective in extracting the true signal. 
1 5 The measurement of oxygen saturation in the blood of a patient is a common 

physiological measurement the accuracy of which may be compromised by the presence 
of noise. Knowledge of blood oxygen saturation can be critical during surgery. There 
are means of obtaining blood oxygen saturation by invasive techniques, such as 
extracting and testing blood removed from a patient using a co-oximeter. But, such 
invasive means are typically time consuming r expensive, and uncomfortable for the 
patient Fortunately, non-invasive measurements of blood oxygen saturation can be 
made using known properties of energy attenuation as a selected form of energy passes 
through a bodily medium. Such non-invasive measurements are performed routinely 
with a pulse oximeter. 

The basic idea behind such energy attenuation measurements is as follows. 
Radiant energy is directed toward a bodily medium, where the medium is derived from 
or contained within a patient, and the amplitude of the energy transmitted through or 
reflected from the medium is then measured. The amount of attenuation of the incident 
energy caused by the medium is strongly dependent on the thickness and composition 
of the medium through which the energy must pass, as well as the specific form of 
energy selected. Information about a physiological system can be derived from data 
taken from the attenuated signal of the incident energy transmitted or reflected. 
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However, the accuracy of such information is reduced where the measured signal 
includes noise. Furthermore, non-invasive measurements often do not afford the 
opportunity to selectively observe the interference causing the noise signal component, 

making it difficult to remove. 

A pulse oximeter is one example of a physiological monitoring system which is 
based upon the measurement of energy attenuated by biological tissues and substances. 
More specifically, a pulse oximeter measures the variable absorption caused by arterial 
blood volume changes. Pulse oximeterettansnut electoomagnetic energy at two 
different wavelengths, typically at 660 nm (red) and 940 nm (infrared, hereinafter IR) 
into the tissue and measure the attenuation of the energy as a function of time. The 
output signal of a pulse oximeter is sensitive to the pulsatile portion of the arterial blood 
flow and contains a component which is a waveform representative of the patient's 
arterial pulse. This type of signal, which contains a component related to the patient's 
pulse, is called a plethysmographic waveform or plethysmogram. 

Pulse oximetry measurements typically use a digit, such as a finger, or an ear 
lobe or other element of the body, where blood flows close to the skin as the medium 
through which light energy is transmitted. The finger, for example, is composed of 
various tissues and substances including skin, fat, bone, muscle, blood, etc. The extent 
to which each of these biological tissues and substances attenuate incident 
electromagnetic energy is generally known. However, the effect of motion can cause 
changes in the optical coupling of the sensor (or probe) to the finger, the underlying 
physiology, the local vasculature, optical properties of tissues due to changing optical 
path length as well as combinations and interactions of the all of the above. Thus, 
patient motion may cause erratic energy attenuation. 

A typical pulse oximeter includes a sensor, cabling from the sensor to a 
computer for signal processing and visual display, the computer and visual display 
typically being included in a patient monitor. The sensor typically includes two light 
emitting diodes (LEDs) placed across a finger tip and a photodetector on the side 
opposite the LEDs. Each LED emits a light signal at different frequencies. The 
detector measures both transmitted light signals once they have passed through the 
finger. The signals are routed to a computer for analysis and display of the various 
parameters measured. 
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The underlying physical basis of a pulse oximeter is Beer's law (also referred to 
as Beer-Lambert's or Bouguer's law) which described attenuation of monochromatic 
light traveling through a uniform medium which absorbs light with the equation: 

/ = 7 £ -dca{X) (1) 

u transmitted 2 incidenT 9 v *' 

where Itransniaed is the intensity of the light transmitted through the uniform medium, 
Jin**/** * s ^ intensity of incident light, d is the distance light is transmitted through the 
uniform medium, c is the concentration of the absorbing substance in the uniform 
medium, expressed in units of mmol L"\ and a(X) is the extinction or absorption 
coefficient of the absorbing substance at wavelength A, expressed in units of L/(mmol 
cm). The properties of Beer's law are valid even if more than one substance absorbs 
light in the medium. Each light absorbing substance contributes its part to the total 
absorbance. 

Each LED emits light at different wavelengths, typically red (centered at about 
660 nm) and IR (centered at about 940 nm) frequency bands. The intensity of light 
transmitted through tissue, I tmsmmedi is different for each wavelength of light emitted by 
the LEDs. Oxyhemoglobin (oxygenated blood) tends to absorb IR light, whereas 
deoxyhemoglobin (deoxygenated blood) tends to absorb red light Thus, the absorption 
of IR light relative to the red light increases with oxyhemoglobin. The ratio of the 
absorption coefficients can be used to determine the oxygen saturation of the blood. 

To estimate blood oxygen saturation, SpO^ a two-solute concentration is 
assumed. A measure of functional blood oxygen saturation level, Sp0 29 can be defined 
as: 

Sp0 2 = °* , (2) 

where c 0 represents oxyhemoglobin solute concentration, and c r represents reduced or 
deoxyhemoglobin solute concentration. 

Noise signal components in a measured pulse oximetry light signal can originate 
from both AC and DC sources. DC noise signal components may be caused by 
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transmission of electromagnetic energy through tissues of relatively constant thickness 
within the body, e.g., bone, muscle, skin, blood, etc. Such DC noise signal components 
may be easily removed with conventional filtering techniques. AC noise signal 
components may occur when tissues being measured are perturbed and, thus, change in 
5 thickness while a measurement is being made. Such AC noise signal components are 
difficult to remove with conventional filtering techniques. Since most materials in and 
derived from the body are easily compressed, the thickness of such matter changes if 
the patient moves during a non-invasive physiological measurement Thus, patient 
movement can cause the properties of energy attenuation to vary erratically. The erratic 

10 or unpredictable nature of motion artifacts induced by noise signal components is a 
major obstacle in removing them. 

Various approaches to removing motion artifacts from measured physiological 
signals, and particularly for use in pulse oximeters, have been proposed. U.S. Patent 
Nos. 5,482,036, 5,490,505, 5,632,272, 5,685,299, and 5,769,785, all to Diab et al., and 

15 U.S. Patent No. 5,91 9,1 34 to Diab, disclose methods and apparatuses for removing 
motion artifacts using adaptive noise cancellation techniques. The basic proposition 
behind the Diab et al. approach is to first generate a noise reference signal from the two 
measured signals, and then use the noise reference signal as an input to an adaptive 
noise canceler along with either or both of the measured signals to remove the reference 

20 noise signal from the measured signals, thus approximating the actual parametric 
signals of interestrT^ 

measured input signals to generate ia noise reference signal. 

Another approach to noise artifact elimination is disclosed in U.S. Patent No. 
5,588,427 to Tien. Tien uses fractal dimension analysis to determine the complexity of 

25 waveforms in order to determine the proper value of the ratio of true intensities based 
on signal complexity. The Tien approach employs a fractal analyzer to determine 
values for two ratios, a and p, based on the measured time varying intensity of the 
transmitted red and IR light signals including noise, a is defined as the ratio of the time 
varying true intensity of light transmitted from the red LED and the time varying true 

30 intensity of the light transmitted from the IR LED. p is a similar ratio relating the noise 
introduced during the measurement of the light transmitted by the red LED and the 
noise introduced during the measurement of the light transmitted by the IRL.ED. 
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According to Tien, a fractal analyzer then determines values for a and p and provides 
(ct,P) pairs to a statistical analyzer. The statistical analyzer performs analysis of one or 
more (a,p) pairs to determine the best value for a, which is then provided to a look-up 
table. The look-up table provides a value corresponding to the arterial oxygen 
saturation in the patient. While the Tien approach appears to be an innovative use of 
fractal analysis, it also appears to be computationally complex. 

Yet another approach to noise artifact elimination is disclosed in U.S. Patent 
Nos. 5,885,213, 5,713,355, 5,555,882 and 5,368,224, all to Richardson et al. The basic 
proposition behind the Richardson et al. approach is to switch operative frequencies 
periodically based on evaluating the noise level associated with various possible 
frequencies of operation in order to select the frequency of operation that has the lowest 
associated noise level. It would appear that data measured at a noisy frequency, using 
the Richardson et al. approach may be invalid or useless for calculating arterial oxygen 
saturation. Furthermore, Richardson et al. requires a computational overhead to 
constantly monitor which frequency of operation provides the least noise. 

Another approach to noise artifact elimination is disclosed in U.S. Patent No. 
5,853,364 to Baker, Jr et al. The Baker, Jr. et al. approach first calculates the heart rate 
of the patient using an adaptive comb filter, power spectrum and pattern matching. 
Once the heart rate is determined, the oximetry data is adaptively comb filtered so that 
only energy at integer multiples of the heart rate are processed. The comb filtered data 
and the raw oximetry data are filtered using a Kalman filter to adaptively modify 
averaging weights and averaging times to attenuate motion artifact noise. The adaptive 
filtering of the Baker, Jr. et al. approach appears to add significant computational 
complexity to solve the problem of motion artifact rejection. 

Still another approach to noise artifact elimination is disclosed in U.S. Patent 
No. 5,431,170 to Mathews. Mathews couples a conventional pulse oximeter light 
transmitter and receiver with a transducer responsive to movement or vibration of the 
body. The transducer provides an electrical signal varying according to the body 
movements or vibrations, which is relatively independent of the blood or other fluid 
flow pulsations. Mathews then provides means for comparing the light signals 
measured with the transducer output and performing adaptive noise cancellation. An 
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apparent disadvantage of the Mathews approach is the need for a secondary sensor to 
detect motion. 

Thus, a need in the art exists for a method, apparatus and system to eliminate 
motion-induced noise artifacts from light signals, that is relatively simple 
5 computationally, and that does not require more than one sensor. 

DISCLOSURE OF INVENTION 
The present invention includes methods, apparatuses and systems for removing 
noise in physiological measurements caused by motion or other similar artifacts. The 
10 methods, apparatuses and systems of the present invention eliminate noise from light 
signals using a single conventional sensor and are relatively simple computationally. 

In a method embodiment, a segment of data from a measured pulse oximetry 
signal is conventionally filtered, and frequency analyzed for major frequency 
components. The frequency components with the largest power spectral density are 
15 selected for subdividing into subsegments, each comprising an individual heartbeat. 
The subsegments are averaged and then analyzed to determine if the averaged 
subsegment is a valid pulse oximetry signal. Additionally, various quality or 
confidence measures may be used to evaluate the validity of such signal. Valid 
averaged subsegments become outputs for further processing to calculate physiological 
20 parameters such as blood oxygen saturation levels. 

A circuit card embodiment includes a processor and memory for storing a 
computer program capable of executing instructions embodying the above method. 

A system embodiment includes an input device, an output device, a memory 
device and a motion artifact rejection circuit card capable of executing instructions 
25 stored in the memory device implementing the methods described herein. 

Finally, a system embodiment includes an input device, and output device, a 
memory device and a processor, which may be a digital signal processor, capable of 
executing instructions stored in the memory device implementing the methods 
described herein. 

30 
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BRIEF DESCRIPTION OF DRAWINGS 
In the drawings, which illustrate what is currently regarded as the best mode for 
carrying out the invention and in which like reference numerals refer to like parts in 
different views or embodiments: 
5 FIG. 1 is a high-level flowchart of a method embodiment of the invention. 

FIG. 2 is a graph of a segment of data representing transmitted infrared (IR) 
light data signal received by a pulse oximetry sensor suitable for signal processing in 
accordance with the invention. 

FIG. 3 is a graph of the power spectrum of the IR data segment in FIG. 2 in 
1 0 accordance with the invention. 

FIG. 4 illustrates a bandpass filtered IR data segment superimposed on the 
original IR data segment of FIG. 2 in accordance with the invention. 

FIG. 5 is a graph of each pulse subsegment detrended and superimposed upon 
one another in accordance with the second method of computing an average pulse. 
15 FIG. 6 shows the modified average pulse as computed from the detrended pulse 

subsegments shown in FIG. 5. 

FIG. 7 illustrates the acquisition of a new pulse segment of raw IR data and the 
removal of the oldest segment of raw IR data in accordance with the invention. 

FIG. 8 is a block diagram of a motion artifact rejection circuit card configured to 
20 remove noise artifacts from signals representing bodily parameters in accordance with 
the invention. 

FIG. 9 is a block diagram of a pulse oximetry system including a motion artifact 
rejection circuit card capable of removing noise from pulse oximetry data in accordance 
with the invention. 

25 FIG. 10 is a block diagram of a pulse oximetry system including a processor 

device programmed to remove noise from pulse oximetiy data in accordance with the 
invention. 

FIGS. 1 1 and 12 are a detailed flowchart of the method of determining the best 
frequency component, block 120, of FIG. 1. 



30 
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BEST MODES FOR CARRYING OUT THE INVENTION 
The following detailed description discloses methods, apparatuses and systems 
for removing motion artifacts from measured plethysmographic waveforms, 
particularly, but without limitation, those used in pulse oximetry. A system 
5 embodiment of the invention includes pulse oximetry hardware and associated software 
to perform the motion artifact suppression. A method embodiment of this invention 
includes a series of steps which exploits certain characteristics of plethysmographic 
waveforms. The methods, apparatuses and systems described below are suitable for 
light transmitted or reflected through bodily tissues and substances. For convenience, 
1 0 the following detailed description will assume measurement of light which has been 
transmitted through a finger of a human. The terms "signal" and "waveform" are used 

interchangeably herein. 

FIG. 1 is a high-level flowchart of a method embodiment of the invention. The 
method steps include acquiring a segment of raw data 100, either red or IR, analyzing 
1 5 the data segment for dominant frequency components 1 1 0, determining the frequency 
component which represents a valid plethysmographic pulse 120, computing an average 
pulse based on the correct frequency component 130 and repeating for new raw data 
segments 140. In order to calculate blood oxygen concentration, Sp0 2 , the method 
embodiment of me invention may be applied to both red and IR data signals to 

20 eliminate or reduce noise from the data signals prior to calculating Sp0 2 . 

The method of this invention begins with acquiring a segment of data (e.g., five 
or more pulses or approximately ten seconds) measured from a single light source 
transmitted through a finger and detected with a sensor on the opposite side of the 
finger. Acquiring a data segment is illustrated by block 100 of FIG. 1. For 

25 convenience, a 1 0.24 second segment of data will be used to illustrate the method. This 
corresponds to 1 024 data points with a sampling rate of 1 00 data points per second. It 
should be readily apparent that the method of the invention is not limited to data 
segments of this size. FIG. 2 is an example of such a data segment for an IR light 
source. The signal processing steps described herein may be performed on both red and 

30 IR data segments independently and simultaneously. Thus, while the steps of the 
method may be illustrated with data from an IR light signal, the same steps are 
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applicable to data from a red light signal and vice versa. The terms "data segment**, 
"input waveform**, "data signal" and "signal** are used interchangeably herein. 

A segment of data may be received from a sensor that converts transmitted or 
reflected light signals into electrical signals. Once a segment of data from a single 
electrical signal has been acquired, it may be filtered to reduce spectral leakage 
resulting from frequency analysis. There are several window filters which may be 
suitable for such purposes. For example, and not by way of limitation, a Hanning 
window may be used to reduce spectral leakage. It will be readily apparent to one of 
skill in the art of digital signal processing that other window filters and methods of 
filtering data to reduce spectral leakage may be selected for reducing spectral leakage. 
As methods of filtering and various filters are known to one of skill in the art of signal 
processing, they will not be further detailed herein. The filtered data is then frequency 
analyzed to determine the dominant frequency components, see block 110 of FIG. 1. 
FIG. 3 illustrates the power spectrum of the IR data segment of FIG. 1 after filtering. 

Signal processing as described herein is generally performed in the frequency 
domain. The segment of data is converted into the frequency domain by, for example, 
performing the well-known Fast Fourier Transform (FFT) on the data. Other common 
techniques of converting time-domain data to the frequency domain may also be used, 
e.g., classical methods using the FFT such as the periodogram or correlogram, 
autoregressive methods, Prony's method, minimum variance methods, maximum 
likelihood methods. Additionally, time domain data may be converted to the frequency 
domain using transforms such as discrete cosine transform, wavelet transform, discrete 
Hartley transform, and Gabor transform. 

Both transient and periodic noise artifacts can induce peaks in the frequency 
domain that may be larger than the peak caused by the patient's heart rate. The 
frequency peak which actually represents the patient's heart rate must then be 
determined, see block 120 of FIG. 1. One approach to determining the correct 
frequency is to order the frequencies by peak amplitude from largest to smallest, F t to 
F c , where F j through F n are not harmonics of each other, and analyze them one by one 
to find the correct frequency, Le., the patient's heart rate. For purposes of illustration, 
only the frequencies associated with the two largest power spectrum amplitude (peaks), 
F, and F 2 , will be used to explain the signal processing in accordance with the 
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invention. It will be readily apparent that the signal processing describe herein may be 
extended from 2 to n candidate frequencies. For convenience of notation, F, is the 
frequency of the largest amplitude peak, and F 2 is the next largest peak, which is not a 
harmonic of F,. 

5 Where this is not the first analysis cycle, an additional check is made to 

determine if one of the two potential frequencies, F, and F 2 , is similar to a known valid 
frequency, F* of the patient's heart rate as determined during the previous analysis 
cycle. Otherwise, the signal processing proceeds as described in the next paragraph. 
Historical trends in heart rate may be used to select the proper frequency peak. In the 

10 case where none of the candidate frequencies is similar to the previous heart rate 

frequency, e.g., both frequencies F, and F 2 are large amplitude noise frequencies, the 
smaller amplitude frequency of the two potential frequencies is discarded, and the 
previous heart rate frequency is selected as the second potential frequency. Thus, the 
method of the invention prevents the situation where there is no correct frequency to 

15 choose because of multiple large amplitude noise frequencies. 

Once candidate frequencies F, and F 2 have been selected, each is processed 
separately to determine which is more likely to be the frequency of the patient's heart 
rate. Analyzing F, first, the data segment may be optionally filtered with a narrow 
bandstop filter at (or near) F 2 to "notch out" that frequency's influence on the data. 

20 Alternatively, each candidate frequency is analyzed in turn without 'Tillering out" the 
effects.oi^ej0Jher_fre(queM determining the 

beginning and end of each heartbeat pulse within the data segment in order to 
subsegment the data into its consecutive heartbeat pulses. 

The period of rhythmic contraction of the heart by which blood is driven 

25 through the aorta and pulmonary artery is known as systole. Maximum light 
absorbance occurs during the systole of a cardiac cycle and is indicated on a 
plethy sinogram by a low point or systolic valley. Conversely, the period of rhythmic 
relaxation and dilation of the heart cavities occurs during diastole when blood is drawn 
into the heart cavities. Minimum light absorbance occurs during the diastole of a 

30 cardiac cycle and is indicated on a plethysmogram by a high point or diastolic peak. 
While it is theoretically possible to use a first derivative of the data segment to identify 
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transitions between diastole and systole and vice versa, in practice, the data may be so 
noisy that the first derivative provides useless information. 

Subsegmenting each pulse within a data segment according to the invention 
begins with a narrow bandpass filter at (or near) F, with a frequency spread (bandpass 
window width) of approximately ± 0.25 Hertz about the center frequency of F„ see 
block 200 of FIG. 1 1 . To improve discrimination, especially with closely spaced peaks, 
the bandpass filter coefficients may be generated and adjusted as needed so that the 
center frequency is nearly identical to the candidate frequency, F,. The resulting 
bandpass filtered data will resemble a sinusoidal waveform and is used to identify the 
point of diastolic peak for each pulse. Diastolic peaks in the data segment will occur at 
or near peaks in the sinusoidal waveform of the bandpass filtered data. FIG. 4 
illustrates the bandpass filtered data superimposed on the original IR data segment of 
FIG. 2. 

Each pulse subsegment is defined as starting one quarter pulse width before a 
diastolic peak and ending approximately one quarter pulse width after the next diastolic 
peak, see block 210 of FIG. 11. This ensures that each heartbeat pulse is centered 
within its pulse subsegment The DC level or bias is determined by taking the mean of 
the data points in each subsegment Each pulse subsegment may be detrended by 
subtracting the DC level, or pulse subsegment mean, from each data point in the pulse 
subsegment 

After subsegmenting the data segment into individual pulses, it is preferred to 
compute the dispersion of saturation values for the latest pulse, see block 220 of FIG. 
11. Measures of the distribution or dispersion of the oxygen saturation values may be 
calculated using the individual data samples, or combinations thereof with a pulse or 
averaged pulse. The measure of dispersion may be used to determine whether there is 
sufficient noise in the data segment to warrant additional signal processing. For 
example, where the dispersion is very low, then the use of previous data in combination 
with the latest pulse would not be required, see blocks 230 and 250 of FIG. 11. Where 
the dispersion is higher, then the pulses for the time interval being evaluated may be 
combined to create an average pulse, see blocks 230 and 240 of FIG. 1 1 . The 
individual pulse or average pulse may then be evaluated using a number of criteria or 
rules, see below. If the pulse under evaluation is found to be acceptable, the saturation 
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could be calculated a number of ways including using the mean or mode of the 
dispersion of the data samples of the individual pulse or average pulse or at a particular 

time or phase of the pulse. 

Additionally, it is important to minimize the dispersion of the saturation values 
5 because the red and IR signals operate on an alternating basis in the time domain, see 
block 270 of FIG. 12. The preferred embodiment minimizes the dispersion of the 
saturation values by adjusting the phase relationship of the red and IR signals. The 
phase relationship of the red and IR signals may be adjusted, for example and not by 
way of limitation, by changing the delay or filtering the signals. Additionally, the delay 

1 0 required adjust the phase relationship of the red and IR signals, and thus, to minimize 
the dispersion of the saturation values, may be variable. It is believed that the variation 
in this delay may be due in part to system to system differences (e.g., monitor and/or 
sensor differences) and other factors (e.g., drive current changes due to physiological 
and anatomical specifics of the patient's measurement site). Because the optimum 

1 5 phase delay between the signals varies, it is desirable to adjust it on a beat to beat basis. 

An exemplary method of irriiumizing dispersion of saturation values varies the 
delay over a range such as about ± 20 ms in 1 ms steps from the nominal values. The 
range may be swept sequentially (i.e., from - 20 ms to +20 ms) or may be performed 
optimally using analytical or empirical search techniques, for example, and not by way 

20 of limitation, by using the Simplex Method or any other suitable search technique. The 
delayed interpolated. waveform may 

order interpolation methods. Methods of interpolation are known to one of skill in the 
art and will not be further detailed. At each step the saturation is calculated for each 
data sample in the individual or average pulse waveform. A histogram of saturation 
25 values is generated using bin sizes of 0.25 percent A measure of the dispersion is then 
calculated. 

The currently preferred measure of dispersion is to find the percentage of all 
saturation values calculated falling within the tallest bin and its two immediately 
surrounding bins. The delay is then adjusted and the process repeated until all steps of 
30 the range have been tested or a minimum dispersion found using an analytical or 

empirical search algorithm. Dispersion may be considered "very low" where at least 
75% of the saturation values calculated fall within the largest (or tallest) bin and its 
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immediately surrounding bins, Le. 9 75% of the calculated saturation values fall in a 
saturation range of 0.75% (3 bin widths). The method of minimizing dispersion by 
locating an optimal delay as described above is exemplary only. Other methods of 
minimizing dispersion and locating an optimal delay will be readily apparent to one of 
skill in the art and should be considered within the spirit and scope of the invention. 

Any suitable means for computing the "central tendency" of the subsegrnented 
pulses or an "average pulse" may be used according to the invention, see block 240 of 
FIG. 1 1 . Two modified mean averaging methods for computing an average pulse will 
be detailed below. Both averaging methods are modified by eliminating data points 
used for computing the average according to various statistical outlier tests. Any 
suitable statistical outlier test may be used according to the invention and those detailed 
below are exemplary only. The terms "average" and "mean" are used interchangeably 
throughout 

The first method for computing an average pulse includes storing the data 
subsegments in a two-dimensional (2-D) array in computer memory so that each row of 
the 2-D array represents a subsegment A modified mean is applied to each of the 
columns of the 2-D array, respectively. The mean is modified by first determining the 
standard deviation of the equivalent points of the subsegments. The standard deviation 
will allow a statistical determination as to whether a point should be included in the 
calculation of the mean or excluded. The preferred embodiment of this statistical 
outlier test is more pessimistic about points at the beginning and end of the 
subsegments (± 1 standard deviation) and less pessimistic about points at the middle of 
the subsegments (± 3 standard deviations). The threshold profile for including data 
points so defined may be linear from beginning to middle and linear from middle to 
end. Alternatively, the threshold profile may have a smooth curvature with maximum 
at the middle of the subsegments. Once the 2-D array has been averaged using the 
modified mean technique, the DC level associated with the last subsegment is applied 
to the averaged subsegment (trended) in the case where the data was previously 
detrended. 

The second method for computing an average pulse includes calculating the 
mean and standard deviation of each pulse subsegment Applying the statistical outlier 
of this second method, only data points within ± 2 standard deviations are used to 
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recalculate a new mean for each pulse subsegment. Each pulse subsegment is 
detrended by subtracting each new mean from each pulse subsegment. Once each pulse 
subsegment is detrended, an ensemble (or column-wise) average of the subsegments 
(/.*., diastolic peak is averaged with diastolic peak, systolic valley is averaged with 
5 systolic valley, etc.) is calculated. This "ensemble" average may be accomplished by 
storing each detrended pulse in a row of a two-dimensional (2-D) array and taking the 
average of the data points in each column. FIG. 5 illustrates twelve individual pulse 
subsegments detrended and superimposed upon one another. FIG. 6 shows the 
modified average pulse as computed from the detrended pulse subsegments shown in 
10 FIG. 5. 

Once the modified average pulse is determined for each candidate frequency, 
they are compared with each other to determine which one most closely matches a 
plethysmographic pulse. This determination is achieved by evaluating various known 
characteristics of a valid plethysmographic pulse. Such characteristics may take the 

1 5 form of rules relating to pulse timing, slopes, etc. that may be used to evaluate the 

validity of a candidate modified average pulse, see block 260 of FIG. 12. For example, 
one can observe the timing between systolic valley and diastolic peak to ensure that the 
time period falls within normal physiological limits. One can assume that the transition 
from diastole to systole is rapid and the transition from the systole back to diastole is 

20 relatively slow. Exemplary rules used to evaluate a candidate modified average pulse 
(or a most recent pulse if. the dispersion is sufficiently .low) include, but are not limited 
to: 

1. The time elapsed from diastolic peak to next diastolic peak,/,, must obey 
physiological limits, for example, 240 milliseconds < /, < 2 seconds, or heart rate, HR, 

25 measured in beats per minute (bpm), must be within physiological limits, for example, 
30 bpm < HR < 250 bpm. 

2. Similarly, the raw signal level at diastolic peak must be less than or equal to 
a maximum value which may be representative of a scaled voltage level. The 
maximum value may be selected from empirical measurements or other known 

30 characteristics of the measurement system. 

3. The time elapsed between diastolic peak and systolic valley, f 2 , must obey the 
following inequality: about 50 milliseconds <t 2 < about 400 milliseconds. 
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4. The time elapsed between diastolic peak and systolic valley must be less than 
the time from systolic valley to the next diastolic peak. In other words, the slope of the 
data signal from diastolic peak to systolic valley is steeper than from systolic valley to 
the next diastolic peak. Or, the absolute value of the slope of the data signal from 
diastolic peak to systolic valley should greater than the absolute value of the slope of 
the data signal from systolic valley to the next diastolic peak. 

5. The width of the candidate pulse cannot vaiy outside a given percentage of 
the average pulse width of each of the pulses within the data segment For example the 
candidate pulse width must be within ± 25% of the average pulse width for the pulses 
identified within the 10.24 seconds of data comprising the data segment 

6. The slope of the data signal from diastolic peak to systolic valley must be 
negative (descending) and the slope of the data signal from systolic valley to the next 
diastolic peak must be positive (ascending). 

The above rules are exemplary only, and not intended to be limiting. 
Furthermore, the various rules selected to evaluate pulse shape may be assigned various 
weights to emphasize relative importance. The above rules may been assigned binary 
values depending on whether the rule is satisfied or not. Alternatively, the rules may be 
assigned weights depending on relative importance. In any case, the above rules may 
be applied to a candidate average pulse and are summed to obtain a quality number. 
The more rules that are satisfied, the higher the quality measure. The resulting quality 
measures for each candidate "frequency may be stati stically analyzed using trend 
analysis and a histogram to determine past performance at these frequencies. Such 
statistical analysis is within the knowledge of one of skill in the art, and thus, will not 
be further detailed herein. 

The above processing is performed for all candidate frequencies, see block 280 
of FIG. 12. The modified average pulse with the highest quality measure is selected, 
see block 290 of FIG. 12. The modified average pulse is then scaled to place each 
diastolic peak at the same signal level to allow the pulses to be appended to one another 
without discontinuity. Once the modified average pulse has been scaled it may be made 
available for further processing, for example, in a blood oxygen saturation calculation, 
see block 130 of FIG. 1. 

The above methods may be integrated into apparatuses and/or systems for 
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calculating blood oxygen saturation. In such pulse oximetry apparatuses and systems, 
the same signal processing as described above is performed on both a red and an JR 
datasignal. The scaled and modified average pulses for red and IR signals may be pre- 
processed and then made available to an oxygen saturation algorithm for calculation of 
arterial oxygen saturation. The above methods may be repeated once another full 
heartbeat pulse of data is collected for both the red and IR signals, see block 1 40 of 
FIG 1 The new pulse of data is added to the red and IR segments, respectively, and 
the oldest pulse of data is removed. FIG. 7 illustrates the acquisition of a new pulse 
segment of raw IR data and the removal of the oldest segment of raw IR data. 

Referring to FIG. 8, the apparatus embodiment of this invention is a motion 
artifact rejection circuit card 10 with an I/O device ll,aprocessor 12 and memory 14 
for storingacomputer programmed algorithm for motion artifact rejection as desenbed 
in the above methods. Processor 12 may be a digital signal processor. I/O device 1 1 
may be ^ circuitry that allows communication to and from external circuitry, for 
example, and not by way of limitation, bus interface circuitry. I/O device 11 may 
include a circuit card edge connector for plugging into a pulse oximetry monitor 
system. Memory 14 may be any solid state electronic memory suitable for storing 
digital data including, for example, computer code and measurement data. 

Referring to FIG. 9, the motion artifact rejection circuit card 10 of FIG. 8 may 
be part of a complete pulse oximetry system 1 6 for eliminating motion-induced noxse 
artifacts in electrical signals (as described in the method embodiments above) and 
calculating and displaying physiological parameters. The pulse oximetry system 16 
also includes an input device 18 and an output device 20. Input device 18 may be a 
pulse oximeter sensor with red and IR LED light sources and a photodetector to convert 
transmitted or reflected light into an electrical signal. Output device 20 may be a 
display device such as a cathode ray tube device, liquid costal display, active matrix 
display or any other suitable display device known to one of skill in the art. 
Alternatively, output device 20 may be a printer for producing a permanent or written 
record such as a laser printer, ink jet printer, thermal printer, dot matrix printer or any 
other suitable printer known to one of skill in the art The pulse oximetry system 16 
may be any pulse oximeter which uses the principles of operation as described above. 
A particular pulse oximeter suitable for use with the circuit card embodiment as 
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described above is the Novametrix Medical Systems, Inc., Model 520A, Pulse 
Oximeter. 

Referring to FIG. 10, a block diagram of a pulse oximetry system 22 including a 
processor device 12, an input device 1 8, an output device 20 and a storage device 24, is 
shown. Input device 1 8 may be a pulse oximeter sensor with red and IR LED light 
sources and a photodetector to convert transmitted or reflected light into an electrical 
signal. Output device 20 may be a display device such as a cathode ray tube device, 
liquid crystal display, active matrix display or any other suitable display device known 
to one of skill in the art Alternatively, output device 20 may be a printer for producing 
a permanent or written record such as a laser printer, ink jet printer, thermal printer, dot 
matrix printer or any other suitable printer known to one of skill in the art Storage 
device 24 may be a disk drive, or any kind of solid state electronic memory device 
suitable for storing digital data including, for example, computer code and measurement 
data. 

Although this invention has been described with reference to particular 
embodiments, the invention is not limited to these described embodiments. Rather, it 
should be understood that the embodiments described herein are merely exemplary and 
that a person skilled in the art may make many variations and modifications without 
departing from the spirit and scope of the invention. All such variations and 
modifications are intended to be included within the scope of the invention as defined 
in the appended claims. 
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CLAIMS 

What is claimed is : 

1. A method of removing motion-induced noise artifacts from a single 
electrical signal representative of a pulse oximetry light signal, comprising: 

5 receiving a segment of raw data spanning a plurality of heartbeats from said single 
electrical signal; 

analyzing said segment of raw data for candidate frequencies, one of which said 

candidate frequencies may be representative of a valid plethysmography pulse; 
analyzing each of said candidate frequencies to determine a best frequency; 
1 0 outputting an average pulse signal computed from said segment of raw data and said 
best frequency; and 
repeating the above steps with a new segment of raw data. 

2. The method of claim 1 , wherein said analyzing said segment of raw data 
1 5 for candidate frequencies includes computing the power spectrum of said segment of 

raw data and identifying n candidate frequencies associated with the n largest amplitude 
peaks in power spectrum which are not harmonics of each other. 



20 



3. The method of claim 2, wherein n equals two. 



4. The method of xlaim3^ further_including comparing said two candidate 

frequencies with a previous known valid frequency and replacing the smaller amplitude 
candidate frequency with said previous known valid frequency if neither of said two 
candidate frequencies is within about ± 0.1 Hertz of said previous known valid 
25 frequency. 

5. The method of claim 2, further including filtering said segment of raw 
data to reduce spectral leakage prior to computing said power spectrum. 
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6. The method of claim 2, wherein said analyzing each of said candidate 
frequencies to determine a best frequency by selecting n candidate frequencies 
associated with the n largest amplitude peaks in power spectrum which are not 
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harmonics of each other includes: 
bandpass filtering said segment of raw data; 

subsegmenting said bandpass filtered segment of raw data into individual heartbeat 
subsegments; 

computing an average heartbeat pulse from said individual heartbeat subsegments; 
evaluating said average heartbeat pulse against known plethysmography pulse 

characteristics to obtain a quality measure for said average heartbeat pulse; 
repeating the above steps for each of said n candidate frequencies; and 
selecting a candidate frequency and average heartbeat pulse with highest quality 

measure to be said best frequency, 

7. The method of claim 6, wherein said bandpass filter is a finite impulse 
response (FIR) filter. 

8. The method of claim 7, wherein the coefficients of said FIR filter are 
adjusted to place the center frequency of said FIR filter at said candidate frequency. 

9. The method of claim 6, wherein said subsegmenting said bandpass 
filtered segment of raw data into individual heartbeat subsegments comprises: 
beginning each heartbeat subsegment at one quarter pulse width before a diastolic peak; 
ending each heartbeat subsegment at one quarter pulse width after the next diastolic 

peak, wherein each heartbeat subsegment overlaps with both a preceding 
heartbeat subsegment and a succeeding heartbeat subsegment, and each 
heartbeat pulse is centered within each heartbeat subsegment 

10. The method of claim 6, wherein said computing an average heartbeat 
pulse from said individual heartbeat subsegments comprises: 

calculating an average for each individual heartbeat subsegment; 
aligning corresponding individual heartbeat subsegment points from beginning to end 
for each of said individual heartbeat subsegments; and 
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calculating a modified average pulse by eliminating all data points outside x standard 
deviations from said average for each individual heartbeat subsegment to obtain 
said average heartbeat pulse. 

5 11. The method of claim 10, wherein x equals about two. 

12. The method of claim 6, further comprising minimizing dispersion of 
saturation values, wherein said minimizing dispersion of saturation values includes 
adjusting the phase relationship of red and IR pulses and calculating a minimum 

10 dispersion of saturation values. 

13. The method of claim 12, wherein adjusting the phase relationship of red 
and IR pulses includes calculating saturation values for a plurality of delays between 
said red and IR signals, wherein said plurality of delays are sequentially swept from 

15 about - 20 ms to about + 20 ms in 1 ms increments. 

14. The method of claim 1 2, wherein minimizing dispersion of saturation 
values includes: 

generating a histogram of saturation values with bin sizes of about y percent calculated 
20 with said delays; 

Jocatmg4i.peakJjiii-and-itU a™ 1 calculating a measure 

of dispersion based on the percentage of saturation values which fall within said 
peak bin and it immediately surrounding z bins; and 
selecting an optimum delay corresponding to measure of dispersion with highest 
25 percentage. 

15. The method of claim 14, wherein y equals 0.25. 

16. The method of claim 14, wherein z equals two. 
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17. A method of removing motion-induced noise artifacts from measured 
pulse oximetry signals prior to calculating functional pulsatile blood oxygen 
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concentration, Sp0 2 , comprising: 

acquiring a red data segment spanning a plurality of heartbeats from oldest to newest; 
frequency analyzing said red data segment for candidate red frequency components, one 

of which may be representative of a valid plethysmography pulse; 
analyzing each of said candidate red frequency components to determine a best red 

frequency component; 
outputting an average red pulse signal computed from said red data segment and said 

best red frequency component; 
repeating the above steps with said new red data segment; and 

substantially simultaneously performing the above steps with respect to an infrared data 
segment. 

18. A method of removing motion-induced noise artifacts from a single 
electrical signal representative of a pulse oximetry light signal comprising: 

(a) acquiring a segment of raw data spanning a plurality of heartbeats from said single 

electrical signal; 

(b) frequency analyzing said segment of raw data to determine two candidate 

frequencies with greatest amplitude, F, and F 2 , which are not harmonics of each 
other; 

(c) removing F 2 from said segment of raw data to obtain F 2 filtered data; 

■(d) subsegmenting said F 2 filtered data into a plurality of subsegments each containing a 
single heartbeat pulse; 

(e) calculating a first average pulse signal based on averaging said plurality of F 2 

filtered subsegments; 

(f) removing F, from said segment of raw data to obtain F, filtered data; 

(g) subsegmenting said F, filtered data into a plurality of subsegments each containing a 

single heartbeat pulse; 

(h) calculating a second average pulse signal based on averaging said plurality of F, 

filtered subsegments; 

(i) determining which of said first average pulse signal and said second average pulse 

signal represents a valid plethysmography pulse; 
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(j) outputting said first average pulse signal or said second average pulse signal which 

represents a valid plethysmographic pulse; and 
(k) repeating steps (b) through (j) with a new segment of raw data. 

5 19. The method of claim 1 8, wherein after step (b), frequencies F, and F 2 are 

compared with a valid frequency, F* from a previous signal processing cycle and F 0 is 
substituted for F 2 if neither F, or F 2 are ± 0.1 Hertz of Fo- 

20. The method of claim 18, wherein removing F 2 from said segment of raw 
10 data is accomplished by filtering said segment of raw data with a bandstop filter of a 

width about 0.5 Hertz centered at F 2 . 

21 . A method of removing motion-induced noise artifacts from a measured 
pulse oximetry signal comprising: 

15 acquiring a segment of pulse oximetry data from the measured pulse oximetry signal; 
filtering said segment of pulse oximetry data; 

performing frequency analysis on said filtered segment of pulse oximetry data to 

determine frequency components of said filtered segment of pulse oximetry 
data; 

20 selecting a frequency component with a largest power as an initial candidate heart rate 
response; 

subsegmenting said filtered segment of pulse oximetry data into a plurality of 

subsegments, wherein each subsegment represents an individual heartbeat; 
averaging data values for each subsegment at identical points along said length of each 
25 subsegment to form an average subsegment; 

determining whether said average subsegment represents a valid pulse oximetry signal; 

outputting valid average subsegments; 

updating said segment of pulse oximetry data; and 

repeating the above steps. 



30 



22. The method of claim 21, further comprising measuring dispersion of 
saturation of latest pulse subsegment to determine whether said averaging step may be 
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skipped. 

23. The method of claim 21, further comprising minimizing dispersion of 
saturation values prior to said outputting valid average subsegments. 

24. A circuit card for use in a pulse oximetry system to remove motion- 
induced noise artifacts from a measured pulse oximetry signal, said circuit card 
comprising: 

a circuit board for mounting electronic circuitry and interfacing with the pulse oximetry 
system; 

a processor mounted on said circuit board for processing at least one input signal 

according to instructions; and 
memory for storing a computer program, wherein said memory is operably coupled to 

said processor, and wherein said computer program includes instructions for 

implementing a method of removing motion artifacts from said measured pulse 

oximetry signal, said method comprising: 

receiving a segment of raw data spanning a plurality of heartbeats from said 

single electrical signal; 
analyzing said segment of raw data for candidate frequencies one of which said 

candidate frequencies may be representative of a valid plethysmography 

pulse; 

analyzing each of said candidate frequencies to determine a best frequency 

based on highest quality measure; 
outputting an average pulse signal computed from said segment of raw data and 

said best frequency; and 
repeating the above steps with a new segment of raw data. 

25. The circuit card of claim 24, wherein said processor is a digital signal 
processor. 



26. A pulse oximetry system for removing motion-induced noise artifacts 
from measured pulse oximetry signals comprising an input device, an output device, 
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and a motion artifact circuit card, wherein said motion artifact circuit card comprises: 
a circuit board configured for coupling to said pulse oximetry system and 

communicating with other electronics within said pulse oximetry system; 
a processor mounted on said circuit board for processing at least one input signal 

according to instructions; and 
memory operably coupled to said processor for storing a computer program, wherein 

said computer program includes instructions for implementing a method of 

removing motion artifacts from said measured pulse oximetry signal, wherein 

said method comprises: 

receiving a segment of raw data spanning a plurality of heartbeats from said 

single electrical signal; 
analyzing said segment of raw data for candidate frequencies one of which said 

candidate frequencies may be representative of a valid plethysmography 

pulse; 

analyzing each of said candidate frequencies to determine a best frequency 

based on highest quality measure; 
outputting an average pulse signal computed from said segment of raw data and 

said best frequency; and 
repeating the above steps with a new segment of raw data. 

27. The pulse-oximetiy -system of claim 26,. wherein_said_processor is a 

digital signal processor. 

28. A pulse oximetry system for removing motion-induced noise artifacts 
from measured pulse oximetry signals comprising an input device, an output device, 
and a motion artifact circuitry, wherein said motion artifact circuitry includes: 
a processor for processing at least one input signal according to instructions; and 
memory operably coupled to said processor for storing a computer program, wherein 

said computer program includes instructions for implementing a method of 

removing motion artifacts from said measured pulse oximetry signal, wherein 

said method comprises: 

receiving a segment of raw data spanning a plurality of heartbeats from said 
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single electrical signal; 
analyzing said segment of raw data for candidate frequencies one of which said 
candidate frequencies may be representative of a valid plethysmographic 
pulse; 

analyzing each of said candidate frequencies to determine a best frequency 

based on highest quality measure; 
outputting an average pulse signal computed from said segment of raw data and 

said best frequency; and 
repeating the above steps with a new segment of raw data. 

29. The pulse oximetry system of claim 28, wherein said processor is a 
digital signal processor. 
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FIG. 2 

Raw Infrared Data Segment 
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FIG, 3 

Power Spectrum of Data from FIG. 1 
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FIG. 4 

Bandpass Filtered and Raw IR Data 
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FIG. 5 
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FIG. 6 
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METHOD, APPARATUS AND SYSTEM FOR REMOVING MOTION 
ARTIFACTS FROM MEASUREMENTS OF BODILY PARAMETERS 

5 PRIORITY CLAIM 

This application claims the benefit of the filing date of United States Provisional 
Patent Application Serial Number 60/104,422, filed October 15, 1998, for "METHOD 
FOR REMOVING MOTION ARTIFACTS FROM DEVICES FOR SENSING 
BODILY PARAMETERS AND APPARATUS FOR EFFECTING SAME". 

10 

TECHNICAL FIELD 
This invention relates to the field of signal processing. More particularly, this 
invention relates to processing measured signals to remove unwanted signal 
components caused by noise and especially noise caused by motion artifacts. 

15 

BACKGROUND 

The measurement of physiological signals can often be difficult because the 
underlying physiological processes may generate very low level signals. Furthermore, 
interfering noise is inherent in the body and the interface between the body and sensors 

20 of physiological processes. Examples of physiological measurements include: 

measurement of electrocardiogram (ECG) signals based on the electrical depolarization 
of the heart muscle, blood pressure, blood oxygen saturation, partial pressure of CO a , 
heart rate, respiration rate, and depth of anesthesia. ECG signals, for example, are 
typically detected by surface electrodes mounted on the chest of a patient ECG signals 

25 are weak at the signal source (/. e. , the heart) and are even weaker at the surface of the 
chest Furthermore, electrical interference from the activity of other muscles {e.g„ 
noise caused by patient breathing, general movement, etc.) causes additional 
interference with physiological signals such as an ECG. Thus, considerable care must 
be taken in the design and use of physiological processors to enhance the quality of the 

30 true signal and reduce the effects of interfering noise signals. 

It is convenient to characterize a measured signal as being a composite signal 
composed of a true signal component and a noise signal component. The terms 
"measured signal" and "composite signal" will be used interchangeably hereinafter. 
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Signal processors are frequently used to remove noise signal components from a 
composite measured signal in order to obtain a signal which closely, if not identically, 
represents the true signal. Conventional filtering techniques such as low pass, band 
pass, and high pass filtering can be used to remove noise signal components from the 
measured composite signal where the noise signal component occupies a frequency 
range outside the true signal component More sophisticated techniques for 
conventional noise filtering include multiple notch filters, which are suitable for use 
where the noise signal component exists at multiple, distinct frequencies, all outside the 
true signal frequency band. 

However, it is often the case that the frequency spectrum of the true and noise 
signal components overlap and that the statistical properties of both signal components 
change with time. More importantly, there are many cases where little is known about 
the noise signal component In such cases, conventional filtering techniques are 
ineffective in extracting the true signal. 

The measurement of oxygen saturation in the blood of a patient is a common 
physiological measurement the accuracy of which may be compromised by the presence 
of noise. Knowledge of blood oxygen saturation can be critical during surgery. There 
are means of obtaining blood oxygen saturation by invasive techniques, such as 
extracting and testing blood removed from a patient using a co-oximeter. But, such 
invasive means are typically time consuming, expensive, and uncomfortable for the 
patient Fortunately, non-invasive measurements of blood oxygen saturation can be 
made using known properties of eneigy attenuation as a selected form of energy passes 
through a bodily medium. Such non-invasive measurements are performed routinely 
with a pulse oximeter. 

The basic idea behind such energy attenuation measurements is as follows. 
Radiant energy is directed toward a bodily medium, where the medium is derived from 
or contained within a patient, and the amplitude of the energy transmitted through or 
reflected from the medium is then measured. The amount of attenuation of the incident 
energy caused by the medium is strongly dependent on the thickness and composition 
of the medium through which the eneigy must pass, as well as the specific form of 
energy selected. Information about a physiological system can be derived from data 
taken from the attenuated signal of the incident energy transmitted or reflected. 
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However, the accuracy of such information is reduced where the measured signal 
includes noise. Furthermore, non-invasive measurements often do not afford the 
opportunity to selectively observe the interference causing the noise signal component, 
making it difficult to remove. 
5 A pulse oximeter is one example of a physiological monitoring system which is 

based upon the measurement of energy attenuated by biological tissues and substances. 
More specifically, a pulse oximeter measures the variable absorption caused by arterial 
blood volume changes. Pulse oximeters transmit electromagnetic energy at two 
different wavelengths, typically at 660 nm (red) and 940 nm (infrared, hereinafter IR) 

1 0 into the tissue and measure the attenuation of the energy as a function of time. The 

output signal of a pulse oximeter is sensitive to the pulsatile portion of the arterial blood 
flow and contains a component which is a waveform representative of the patient's 
arterial pulse. This type of signal, which contains a component related to the patient's 
pulse, is called a plethysmography waveform or plethysmogram. 

1 5 Pulse oximetry measurements typically use a digit, such as a finger, or an ear 

lobe or other element of the body, where blood flows close to the skin as the medium 
through which light energy is transmitted. The finger, for example, is composed of 
various tissues and substances including skin, fat, bone, muscle, blood, etc. The extent 
to which each of these biological tissues and substances attenuate incident 

20 electromagnetic energy is generally known. However, the effect of motion can cause 

Tz^T^changejtra^ ~ 
physiology, the local vasculature, optical properties of tissues due to changing optical 
path length as well as combinations and interactions of the all of the above. Thus, 
patient motion may cause erratic energy attenuation. 

25 A typical pulse oximeter includes a sensor, cabling from the sensor to a 

computer for signal processing and visual display, the computer and visual display 
typically being included in a patient monitor. The sensor typically includes two light 
emitting diodes (LEDs) placed across a finger tip and a photodetector on the side 
opposite the LEDs. Each LED emits a light signal at different frequencies. The 

30 detector measures both transmitted light signals once they have passed through the 
finger. The signals are routed to a computer for analysis and display of the various 
parameters measured. 
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The underlying physical basis of a pulse oximeter is Beer's law (also referred to 
as Beer-Lambert's or Bouguer's law) which described attenuation of monochromatic 
light traveling through a uniform medium which absorbs light with the equation: 



where I tranS mitted * s Ae intensity of the light transmitted through the uniform medium, 
Imcidem IS the intensity of incident light, d is the distance light is transmitted through the 
uniform medium, c is the concentration of the absorbing substance in the uniform 
medium, expressed in units of mmol L" 1 , and a(X) is the extinction or absorption 
coefficient of the absorbing substance at wavelength X, expressed in units of L/(mmol 
cm). The properties of Beer's law are valid even if more than one substance absorbs 
light in the medium. Each light absorbing substance contributes its part to the total 
absorbance. 

Each LED emits light at different wavelengths, typically red (centered at about 
660 nra) and IR (centered at about 940 nm) frequency bands. The intensity of light 
transmitted through tissue, I tmsmaud> is different for each wavelength of light emitted by 
the LEDs. Oxyhemoglobin (oxygenated blood) tends to absorb IR light, whereas 
deoxyhemoglobin (deoxygenated blood) tends to absorb red light Thus, the absorption 
of IR light relative to the red light increases with oxyhemoglobin. The ratio of the 
absorption coefficients can be used to determine the oxygen saturation of the blood. 

To estimate blood oxygen saturation, SpO^ a two-solute concentration is 
assumed. A measure of functional blood oxygen saturation level, SpO^ can be defined 
as: 



SpO z = — C -2— , (2) 

where c 0 represents oxyhemoglobin solute concentration, and c r represents reduced or 
deoxyhemoglobin solute concentration. 

Noise signal components in a measured pulse oximetry light signal can originate 
from both AC and DC sources. DC noise signal components may be caused by 
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transmission of electromagnetic energy through tissues of relatively constant thickness 
within the body, e.g., bone, muscle, skin, blood, etc. Such DC noise signal components 
may be easily removed with conventional filtering techniques. AC noise signal 
components may occur when tissues being measured are perturbed and, thus, change in 
5 thickness while a measurement is being made. Such AC noise signal components are 
difficult to remove with conventional filtering techniques. Since most materials in and 
derived from the body are easily compressed, the thickness of such matter changes if 
the patient moves during a non-invasive physiological measurement Thus, patient 
movement can cause the properties of energy attenuation to vary erratically . The erratic 

10 or unpredictable; nature of motion artifacts induced by noise signal components is a 
major obstacle in removing them. 

Various approaches to removing motion artifacts from measured physiological 
signals, and particularly for use in pulse oximeters, have been proposed. U.S. Patent 
Nos. 5,482,036, 5,490,505, 5,632,272, 5,685,299, and 5,769,785, all to Diab et al., and 

15 U.S. Patent No. 5,919,134 to Diab, disclose methods and apparatuses for removing 
motion artifacts using adaptive noise cancellation techniques. The basic proposition 
behind the Diab et al. approach is to first generate a noise reference signal from the two 
measured signals, and then use the noise reference signal as an input to an adaptive 
noise canceler along with either or both of the measured signals to remove the reference 

20 noise signal from the measured signals, thus approximating the actual parametric 
signals of interest The Diab et al. approach appears to require the use of both 
measured input signals to generate a noise reference signal. 

Another approach to noise artifact elimination is disclosed in U.S. Patent No. 
5,588,427 to Tien. Tien uses fractal dimension analysis to determine the complexity of 

25 waveforms in order to determine the proper value of the ratio of true intensities based 
on signal complexity. The Tien approach employs a fractal analyzer to determine 
values for two ratios, a and p, based on the measured time varying intensity of the 
transmitted red and IR light signals including noise, a is defined as the ratio of the time 
varying true intensity of light transmitted from the red LED and the time varying true 

30 intensity of the light transmitted from the IR LED. P is a similar ratio relating the noise 
introduced during the measurement of the light transmitted by the red LED and the 
noise introduced during the measurement of the light transmitted by the IR LED. 
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According to Tien, a fractal analyzer then determines values for a and P and provides 
(a,p) pairs to a statistical analyzer. The statistical analyzer performs analysis of one or 
more (ct,P) pairs to determine the best value for a, which is then provided to a look-up 
table. The look-up table provides a value corresponding to the arterial oxygen 
saturation in the patient While the Tien approach appears to be an innovative use of 
fractal analysis, it also appears to be computationally complex. 

Yet another approach to noise artifact elimination is disclosed in U.S. Patent 
Nos. 5,885,213, 5,713,355, 5,555,882 and 5,368,224, all to Richardson et al. The basic 
proposition behind the Richardson et al. approach is to switch operative frequencies 
periodically based on evaluating the noise level associated with various possible 
frequencies of operation in order to select the frequency of operation that has the lowest 
associated noise level. It would appear that data measured at a noisy frequency, using 
the Richardson et al. approach may be invalid or useless for calculating arterial oxygen 
saturation. Furthermore, Richardson et al. requires a computational overhead to 
constantly monitor which frequency of operation provides the least noise. 

Another approach to noise artifact elimination is disclosed in U.S. Patent No. 
5,853,364 to Baker, Jr et al. The Baker, Jr. et al. approach first calculates the heart rate 
of the patient using an adaptive comb filter, power spectrum and pattern matching. 
Once the heart rate is determined, the oximetry data is adapti vely comb filtered so that 
only energy at integer multiples of the heart rate are processed. The comb filtered data 
and the raw oximetry data are filtered using a Kalman filter to adaptively modify 
averaging weights and averaging times to attenuate motion artifact noise. The adaptive 
filtering of the Baker, Jr. et al. approach appears to add significant computational 
complexity to solve the problem of motion artifact rejection. 

Still another approach to noise artifact elimination is disclosed in U.S. Patent 
No. 5,43 1,170 to Mathews. Mathews couples a conventional pulse oximeter light 
transmitter and receiver with a transducer responsive to movement or vibration of the 
body. The transducer provides an electrical signal varying according to the body 
movements or vibrations, which is relatively independent of the blood or other fluid 
flow pulsations. Mathews then provides means for comparing the light signals 
measured with the transducer output and performing adaptive noise cancellation. An 



WO 00/22408 PCT/US99/23803 

-7- 

apparent disadvantage of the Mathews approach is the need for a secondaiy sensor to 
detect motion. 

Thus, a need in the art exists for a method, apparatus and system to eliminate 
motion-induced noise artifacts from light signals, that is relatively simple 
5 computationally, and that does not require more than one sensor. 

DISCLOSURE OF INVENTION 
The present invention includes methods, apparatuses and systems for removing 
noise in physiological measurements caused by motion or other similar artifacts. The 
10 methods, apparatuses and systems of the present invention eliminate noise from light 
signals using a single conventional sensor and are relatively simple computationally. 

In a method embodiment, a segment of data from a measured pulse oximetry 
signal is conventionally filtered, and frequency analyzed for major frequency 
components. The frequency components with the largest power spectral density are 
15 selected for subdividing into subsegments, each comprising an individual heartbeat 
The subsegments are averaged and then analyzed to determine if the averaged 
subsegment is a valid pulse oximetry signal. Additionally, various quality or 
confidence measures may be used to evaluate the validity of such signal. Valid 
averaged subsegments become outputs for further processing to calculate physiological 
20 parameters such as blood oxygen saturation levels. 

A circuit card embodiment jncl^^ a for storing a 

computer program capable of executing instructions embodying the above method. 

A system embodiment includes an input device, an output device, a memory 
device and a motion artifact rejection circuit card capable of executing instructions 
25 stored in the memory device implementing the methods described herein. 

Finally, a system embodiment includes an input device, and output device, a 
memory device and a processor, which may be a digital signal processor, capable of 
executing instructions stored in the memory device implementing the methods 
described herein. 
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BRIEF DESCRIPTION OF DRAWINGS 
In the drawings, which illustrate what is currently regarded as the best mode for 
carrying out the invention and in which like reference numerals refer to like parts in 
different views or embodiments: 

FIG. 1 is a high-level flowchart of a method embodiment of the invention. 
FIG. 2 is a graph of a segment of data representing transmitted infrared (IR) 
light data signal received by a pulse oximetry sensor suitable for signal processing in 
accordance with the invention. 

FIG. 3 is a graph of the power spectrum of the IR data segment in FIG. 2 in 
accordance with the invention. 

FIG. 4 illustrates a bandpass filtered IR data segment superimposed on the 
original IR data segment of FIG. 2 in accordance with the invention. 

FIG. 5 is a graph of each pulse subsegment detrended and superimposed upon 
one another in accordance with the second method of computing an average pulse. 

FIG. 6 shows the modified average pulse as computed from the detrended pulse 
subsegments shown in FIG. 5. 

FIG. 7 illustrates the acquisition of a new pulse segment of raw IR data and the 
removal of the oldest segment of raw IR data in accordance with the invention. 

FIG. 8 is a block diagram of a motion artifact rejection circuit card configured to 
remove noise artifacts from signals representing bodily parameters in accordance with 
the invention. 

FIG. 9 is a block diagram of a pulse oximetry system including a motion artifact 
rejection circuit card capable of removing noise from pulse oximetry data in accordance 
with the invention. 

FIG. 10 is a block diagram of a pulse oximetry system including a processor 
device programmed to remove noise from pulse oximetry data in accordance with the 
invention. 

FIGS. 1 1 and 12 are a detailed flowchart of the method of determining the best 
frequency component, block 120, of FIG. 1. 
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BEST MODES FOR CARRYING OUT THE INVENTION 
The following detailed description discloses methods, apparatuses and systems 
for removing motion artifacts from measured plethysmographic waveforms, 
particularly, but without limitation, those used in pulse oximetry. A system 
5 embodiment of the invention includes pulse oximetry hardware and associated software 
to perform the motion artifact suppression. A method embodiment of this invention 
includes a series of steps which exploits certain characteristics of plethysmographic 
waveforms. The methods, apparatuses and systems described below are suitable for 
light transmitted or reflected through bodily tissues and substances. For convenience, 

10 the following detailed description will assume measurement of light which has been 
transmitted through a finger of a human. The terms "signal" and "waveform" are used 
interchangeably herein. 

FIG. 1 is a high-level flowchart of a method embodiment of the invention. The 
method steps include acquiring a segment of raw data 100, either red or IR, analyzing 

1 5 the data segment for dominant frequency components 1 1 0, determining the frequency 
component which represents a valid plethysmographic pulse 120, computing an average 
pulse based on the correct frequency component 130 and repeating for new raw data 
segments 140. In order to calculate blood oxygen concentration, Sp0 2 , the method 
embodiment of the invention may be applied to both red and IR data signals to 

20 eliminate or reduce noise from the data signals prior to calculating Sp0 2 . 

The method of this invention begins with acquiring a segment of data {e.g., five 
or more pulses or approximately ten seconds) measured from a single light source 
transmitted through a finger and detected with a sensor on the opposite side of the 
finger. Acquiring a data segment is illustrated by block 100 of FIG. 1. For 

25 convenience, a 1 0.24 second segment of data will be used to illustrate the method. This 
corresponds to 1024 data points with a sampling rate of 1 00 data points per second. It 
should be readily apparent that the method of the invention is not limited to data 
segments of this size. FIG. 2 is an example of such a data segment for an IR light 
source. The signal processing steps described herein may be performed on both red and 

30 IR data segments independently and simultaneously. Thus, while the steps of the 
method may be illustrated with data from an IR light signal, the same steps are 
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applicable to data from a red light signal and vice versa. The terms "data segment", 
"input waveform", "data signal" and "signal" are used interchangeably herein. 

A segment of data may be received from a sensor that converts transmitted or 
reflected light signals into electrical signals. Once a segment of data from a single 
electrical signal has been acquired, it may be filtered to reduce spectral leakage 
resulting from frequency analysis. There are several window filters which may be 
suitable for such purposes. For example, and not by way of limitation, a Harming 
window may be used to reduce spectral leakage. It will be readily apparent to one of 
skill in the art of digital signal processing that other window filters and methods of 
filtering data to reduce spectral leakage may be selected for reducing spectral leakage. 
As methods of filtering and various filters are known to one of skill in the art of signal 
processing, they will not be further detailed herein. The filtered data is then frequency 
analyzed to determine the dominant frequency components, see block 110 of FIG. 1 . 
FIG. 3 illustrates the power spectrum of the IR data segment of FIG. 1 after filtering. 

Signal processing as described herein is generally performed in the frequency 
domain. The segment of data is converted into the frequency domain by, for example, 
performing the well-known Fast Fourier Transform (FFT) on the data. Other common 
techniques of converting time-domain data to the frequency domain may also be used, 
e.g., classical methods using the FFT such as the periodogram or correlogram, 
autoregressive methods, Prony's method, minimum variance methods, maximum 
likelihood methods. Additionally, time domain data may be converted to the frequency 
domain using transforms such as discrete cosine transform, wavelet transform, discrete 
Hartley transform, and Gabor transform. 

Both transient and periodic noise artifacts can induce peaks in the frequency 
domain that may be larger than the peak caused by the patient's heart rate. The 
frequency peak which actually represents the patient's heart rate must then be 
determined, see block 120 of FIG. 1. One approach to determining the correct 
frequency is to order the frequencies by peak amplitude from largest to smallest, F, to 
F n , where F, through F n are not harmonics of each other, and analyze them one by one 
to find the correct frequency, Le. 9 the patient's heart rate. For purposes of illustration, 
only the frequencies associated with the two largest power spectrum amplitude (peaks), 
F, and F 2 , will be used to explain the signal processing in accordance with the 
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invention. It will be readily apparent that the signal processing describe herein may be 
extended from 2 to n candidate frequencies. For convenience of notation, F, is the 
frequency of the largest amplitude peak, and F 2 is the next largest peak, which is not a 
harmonic of F,. 

5 Where this is not the first analysis cycle, an additional check is made to 

determine if one of the two potential frequencies, F, and F 2? is similar to a known valid 
frequency, Fq, of the patient's heart rate as determined during the previous analysis 
cycle. Otherwise, the signal processing proceeds as described in the next paragraph. 
Historical trends in heart rate may be used to select the proper frequency peak. In the 

10 case where none of the candidate frequencies is similar to the previous heart rate 

frequency, e.g., both frequencies F, and F 2 are large amplitude noise frequencies, the 
smaller amplitude frequency of the two potential frequencies is discarded, and the 
previous heart rate frequency is selected as the second potential frequency. Thus, the 
method of the invention prevents the situation where there is no correct frequency to 

15 choose because of multiple large amplitude noise frequencies. 

Once candidate frequencies F, and F 2 have been selected, each is processed 
separately to determine which is more likely to be the frequency of the patient's heart 
rate. Analyzing F, first, the data segment may be optionally filtered with a narrow 
bandstop filter at (or near) F 2 to "notch out" that frequency * s influence on the data. 

20 Alternatively, each candidate frequency is analyzed in turn without "filtering out" the 
effects of the other frequencies. The signal processing continues by determining the 
beginning and end of each heartbeat pulse within the data segment in order to 
subsegment the data into its consecutive heartbeat pulses. 

The period of rhythmic contraction of the heart by which blood is driven 

25 through the aorta and pulmonary artery is known as systole. Maximum light 
absorbance occurs during the systole of a cardiac cycle and is indicated on a 
plethysmogram by a low point or systolic valley. Conversely, the period of rhythmic 
relaxation and dilation of the heart cavities occurs during diastole when blood is drawn 
into the heart cavities. Minimum light absorbance occurs during the diastole of a 

30 cardiac cycle and is indicated on a plethysmogram by a high point or diastolic peak. 
While it is theoretically possible to use a first derivative of the data segment to identify 
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transitions between diastole and systole and vice versa, in practice, the data may be so 
noisy that the first derivative provides useless information. 

Subsegmenting each pulse within a data segment according to the invention 
begins with a narrow bandpass filter at (or near) F ( with a frequency spread (bandpass 
window width) of approximately ± 0.25 Hertz about the center frequency of F„ see 
block 200 of FIG. 1 1. To improve discrimination, especially with closely spaced peaks, 
the bandpass filter coefficients may be generated and adjusted as needed so that the 
center frequency is nearly identical to the candidate frequency, F,. The resulting 
bandpass filtered data will resemble a sinusoidal waveform and is used to identify the 
point of diastolic peak for each pulse. Diastolic peaks in the data segment will occur at 
or near peaks in the sinusoidal waveform of the bandpass filtered data. FIG. 4 
illustrates the bandpass filtered data superimposed on the original IR data segment of 
FIG. 2. 

Each pulse subsegment is defined as starting one quarter pulse width before a 
diastolic peak and ending approximately one quarter pulse width after the next diastolic 
peak, see block 2 1 0 of FIG. 1 1 . This ensures that each heartbeat pulse is centered 
within its pulse subsegment The DC level or bias is determined by taking the mean of 
the data points in each subsegment Each pulse subsegment may be detrended by 
subtracting the DC level, or pulse subsegment mean, from each data point in the pulse 
subsegment 

After subsegmenting the data segment into individual pulses, it is preferred to 
compute the dispersion of saturation values for the latest pulse, see block 220 of FIG. 
1 1 . Measures of the distribution or dispersion of the oxygen saturation values may be 
calculated using the individual data samples, or combinations thereof with a pulse or 
averaged pulse. The measure of dispersion may be used to determine whether there is 
sufficient noise in the data segment to warrant additional signal processing. For 
example, where the dispersion is very low, then the use of previous data in combination 
with the latest pulse would not be required, see blocks 230 and 250 of FIG. 1 1 . Where 
the dispersion is higher, then the pulses for the time interval being evaluated may be 
combined to create an average pulse, see blocks 230 and 240 of FIG. 1 1 . The 
individual pulse or average pulse may then be evaluated using a number of criteria or 
rules, see below. If the pulse under evaluation is found to be acceptable, the saturation 
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could be calculated a number of ways including using the mean or mode of the 
dispersion of the data samples of the individual pulse or average pulse or at a particular 
time or phase of the pulse. 

Additionally, it is important to minimize the dispersion of the saturation values 
5 because the red and IR signals operate on an alternating basis in the time domain, see 
block 270 of FIG. 12. The preferred embodiment minimizes the dispersion of the 
saturation values by adjusting the phase relationship of the red and IR signals. The 
phase relationship of the red and IR signals may be adjusted, for example and not by 
way of limitation, by changing the delay or filtering the signals. Additionally, the delay 

1 0 required adjust the phase relationship of the red and IR signals, and thus, to minimize 
the dispersion of the saturation values, may be variable. It is believed that the variation 
in this delay may be due in part to system to system differences (e.g., monitor and/or 
sensor differences) and other factors (e.g., drive current changes due to physiological 
and anatomical specifics of the patient's measurement site). Because the optimum 

1 5 phase delay between the signals varies, it is desirable to adjust it on a beat to beat basis. 

An exemplary method of minimizing dispersion of saturation values varies the 
delay over a range such as about ± 20 ms in 1 ms steps from the nominal values. The 
range may be swept sequentially (Le. , from - 20 ms to +20 ms) or may be performed 
optimally using analytical or empirical search techniques, for example, and not by way 

20 of limitation, by using the Simplex Method or any other suitable search technique. The 
delayedlntejrp^ may be estimated using linear interpolation or higher 

order interpolation methods. Methods of interpolation are known to one of skill in the 
art and will not be further detailed. At each step the saturation is calculated for each 
data sample in the individual or average pulse waveform. A histogram of saturation 

25 values is generated using bin sizes of 0.25 percent. A measure of the dispersion is then 
calculated. 

The currently preferred measure of dispersion is to find the percentage of all 
saturation values calculated falling within the tallest bin and its two immediately 
surrounding bins. The delay is then adjusted and the process repeated until all steps of 
30 the range have been tested or a minimum dispersion found using an analytical or 

empirical search algorithm. Dispersion may be considered 4t very low" where at least 
75% of the saturation values calculated fall within the largest (or tallest) bin and its 
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immediately surrounding bins, i.e., 75% of the calculated saturation values fall in a 
saturation range of 0.75% (3 bin widths). The method of minimizing dispersion by 
locating an optimal delay as described above is exemplary only. Other methods of 
minimizing dispersion and locating an optimal delay will be readily apparent to one of 
skill in the art and should be considered within the spirit and scope of the invention. 

Any suitable means for computing the "central tendency" of the subsegmented 
pulses or an "average pulse" may be used according to the invention, see block 240 of 
FIG. 1 1 . Two modified mean averaging methods for computing an average pulse will 
be detailed below. Both averaging methods are modified by eliminating data points 
used for computing the average according to various statistical outlier tests. Any 
suitable statistical outlier test may be used according to the invention and those detailed 
below are exemplary only. The terms "average" and "mean" are used interchangeably 
throughout 

The first method for computing an average pulse includes storing the data 
subsegments in a two-dimensional (2-D) array in computer memory so that each row of 
the 2-D array represents a subsegment A modified mean is applied to each of the 
columns of the 2-D array, respectively. The mean is modified by first determining the 
standard deviation of the equivalent points of the subsegments. The standard deviation 
will allow a statistical determination as to whether a point should be included in the 
calculation of the mean or excluded. The preferred embodiment of this statistical 
outlier test is more pessimistic about points at the beginning and end of the 
subsegments (± 1 standard deviation) and less pessimistic about points at the middle of 
the subsegments (± 3 standard deviations). The threshold profile for including data 
points so defined may be linear from beginning to middle and linear from middle to 
end. Alternatively, the threshold profile may have a smooth curvature with maximum 
at the middle of the subsegments. Once the 2-D array has been averaged using the 
modified mean technique, the DC level associated with the last subsegment is applied 
to the averaged subsegment (trended) in the case where the data was previously 
detrended. 

The second method for computing an average pulse includes calculating the 
mean and standard deviation of each pulse subsegment Applying the statistical outlier 
of this second method, only data points within ± 2 standard deviations are used to 
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recalculate a new mean for each pulse subsegment. Each pulse subsegment is 
detrended by subtracting each new mean from each pulse subsegment. Once each pulse 
subsegment is detrended, an ensemble (or column-wise) average of the subsegments 
(j.e., diastolic peak is averaged with diastolic peak, systolic valley is averaged with 
5 systolic valley, etc.) is calculated. This "ensemble" average may be accomplished by 
storing each detrended pulse in a row of a two-dimensional (2-D) array and taking the 
average of the data points in each column. FIG. 5 illustrates twelve individual pulse 
subsegments detrended and superimposed upon one another. FIG. 6 shows the 
modified average pulse as computed from the detrended pulse subsegments shown in 
10 FIG. 5. 

Once the modified average pulse is determined for each candidate frequency, 
they are compared with each other to determine which one most closely matches a 
plethysmographic pulse. This determination is achieved by evaluating various known 
characteristics of a valid plethysmographic pulse. Such characteristics may take the 

1 5 form of rules relating to pulse timing, slopes, etc. that may be used to evaluate the 

validity of a candidate modified average pulse, see block 260 of FIG. 12. For example, 
one can observe the timing between systolic valley and diastolic peak to ensure that the 
time period falls within normal physiological limits. One can assume that the transition 
from diastole to systole is rapid and the transition from the systole back to diastole is 

20 relatively slow. Exemplary rules used to evaluate a candidate modified average pulse 
(or a most recent pulse if the dispersion is sufficiently low) include, but are not limited 
to: 

1 . The time elapsed from diastolic peak to next diastolic peak, f „ must obey 
physiological limits, for example, 240 milliseconds < t x < 2 seconds, or heart rate, HR, 

25 measured in beats per minute (bpm), must be within physiological limits, for example, 
30 bpm < HR < 250 bpm. 

2. Similarly, the raw signal level at diastolic peak must be less than or equal to 
a maximum value which may be representative of a scaled voltage level. The 
maximum value may be selected from empirical measurements or other known 

30 characteristics of the measurement system. 

3. The time elapsed between diastolic peak and systolic valley, r 2 , must obey the 
following inequality: about 50 milliseconds <t 2 < about 400 milliseconds. 
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4. The time elapsed between diastolic peak and systolic valley must be less than 
the time from systolic valley to the next diastolic peak. In other words, the slope of the 
data signal from diastolic peak to systolic valley is steeper than from systolic valley to 
the next diastolic peak. Or, the absolute value of the slope of the data signal from 
diastolic peak to systolic valley should greater than the absolute value of the slope of 
the data signal from systolic valley to the next diastolic peak. 

5. The width of the candidate pulse cannot vary outside a given percentage of 
the average pulse width of each of the pulses within the data segment For example the 
candidate pulse width must be within ± 25% of the average pulse width for the pulses 
identified within the 10.24 seconds of data comprising the data segment 

6. The slope of the data signal from diastolic peak to systolic valley must be 
negative (descending) and the slope of the data signal from systolic valley to the next 
diastolic peak must be positive (ascending). 

The above rules are exemplary only, and not intended to be limiting. 
Furthermore, the various rules selected to evaluate pulse shape may be assigned various 
weights to emphasize relative importance. The above rules may been assigned binary 
values depending on whether the rule is satisfied or not Alternatively, the rules may be 
assigned weights depending on relative importance. In any case, the above rules may 
be applied to a candidate average pulse and are summed to obtain a quality number. 
The more rules that are satisfied, the higher the quality measure. The resulting quality 
measures for each candidate frequency may be statistically analyzed using trend 
analysis and a histogram to determine past performance at these frequencies. Such 
statistical analysis is within the knowledge of one of skill in the art, and thus, will not 
be further detailed herein. 

The above processing is performed for all candidate frequencies, see block 280 
of FIG. 12. The modified average pulse with the highest quality measure is selected, 
see block 290 of FIG. 12. The modified average pulse is then scaled to place each 
diastolic peak at the same signal level to allow the pulses to be appended to one another 
without discontinuity. Once the modified average pulse has been scaled it may be made 
available for further processing, for example, in a blood oxygen saturation calculation, 
see block 130 of FIG. 1. 

The above methods may be integrated into apparatuses and/or systems for 
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calculating blood oxygen saturation. In such pulse oximetry apparatuses and systems, 
the same signal processing as described above is performed on both a red and an 1R 
data signal. The scaled and modified average pulses for red and IR signals may be pre- 
processed and then made available to an oxygen saturation algorithm for calculation of 
5 arterial oxygen saturation. The above methods may be repeated once another full 
heartbeat pulse of data is collected for both the red and IR signals, see block 140 of 
FIG. 1 . The new pulse of data is added to the red and IR segments, respectively, and 
the oldest pulse of data is removed. FIG. 7 illustrates the acquisition of a new pulse 
segment of raw IR data and the removal of the oldest segment of raw IR data. 

1 o Referring to FIG. 8, the apparatus embodiment of this invention is a motion 

artifact rejection circuit card 10 with an I/O device 1 1, a processor 12 and memory 14 
for storing a computer programmed algorithm for motion artifact rejection as described 
in the above methods. Processor 12 may be a digital signal processor. I/O device 1 1 
may be any circuitry that allows communication to and from external circuitry, for 

1 5 example, and not by way of limitation, bus interface circuitry. I/O device 1 1 may 
include a circuit card edge connector for plugging into a pulse oximetry monitor 
system. Memory 14 may be any solid state electronic memory suitable for storing 
digital data including, for example, computer code and measurement data. 

Referring to FIG. 9, the motion artifact rejection circuit card 10 of FIG. 8 may 

20 be part of a complete pulse oximetry system 16 for eliminating motion-induced noise 
artifacts in electrical signals (as described in the method embodiments above) and 
calculating and displaying physiological parameters. The pulse oximetry system 16 
also includes an input device 18 and an output device 20. Input device 18 may be a 
pulse oximeter sensor with red and IR LED light sources and a photodetector to convert 

25 transmitted or reflected light into an electrical signal. Output device 20 may be a 

display device such as a cathode ray tube device, liquid crystal display, active matrix 
display or any other suitable display device known to one of skill in the art 
Alternatively, output device 20 may be a printer for producing a permanent or written 
record such as a laser printer, ink jet printer, thermal printer, dot matrix printer or any 

30 other suitable printer known to one of skill in the art. The pulse oximetry system 16 
may be any pulse oximeter which uses the principles of operation as described above. 
A particular pulse oximeter suitable for use with the circuit card embodiment as 
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described above is the Novametrix Medical Systems, Inc., Model 520A, Pulse 
Oximeter. 

Referring to FIG. 10, a block diagram of a pulse oximetry system 22 including a 
processor device 12, an input device 1 8, an output device 20 and a storage device 24, is 
shown. Input device 1 8 may be a pulse Oximeter sensor with red and IR LED light 
sources and a photodetector to convert transmitted or reflected light into an electrical 
signal. Output device 20 may be a display device such as a cathode ray tube device, 
liquid crystal display, active matrix display or any other suitable display device known 
to one of skill in the art Alternatively, output device 20 may be a printer for producing 
a permanent or written record such as a laser printer, ink jet printer, thermal printer, dot 
matrix printer or any other suitable printer known to one of skill in the art Storage 
device 24 may be a disk drive, or any kind of solid state electronic memory device 
suitable for storing digital data including, for example, computer code and measurement 
data. 

Although this invention has been described with reference to particular 
embodiments, the invention is not limited to these described embodiments. Rather, it 
should be understood that the embodiments described herein are merely exemplary and 
that a person skilled in the art may make many variations and modifications without 
departing from the spirit and scope of the invention. All such variations and 
modifications are intended to be included within the scope of the invention as defined 
in the appended claims. 
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What is claimed is : 

1 . A method of removing motion-induced noise artifacts from a single 
electrical signal representative of a pulse oximetry light signal, comprising: 

5 receiving a segment of raw data spanning a plurality of heartbeats from said single 
electrical signal; 

analyzing said segment of raw data for candidate frequencies, one of which said 

candidate frequencies may be representative of a valid plethysmographic pulse; 

analyzing each of said candidate frequencies to determine a best frequency; 
10 outputting an average pulse signal computed from said segment of raw data and said 
best frequency; and 

repeating the above steps with a new segment of raw data. 

2. The method of claim 1 , wherein said analyzing said segment of raw data 
15 for candidate frequencies includes computing the power spectrum of said segment of 

raw data and identifying n candidate frequencies associated with the n largest amplitude 
peaks in power spectrum which are not harmonics of each other. 



20 



3. The method of claim 2, wherein n equals two. 



4. The method of claim 3, further including comparing said two candidate 
frequencies with a previous known valid frequency and replacing the smaller amplitude 
candidate frequency with said previous known valid frequency if neither of said two 
candidate frequencies is within about ±0.1 Hertz of said previous known valid 

25 frequency. 

5. The method of claim 2, further including filtering said segment of raw 
data to reduce spectral leakage prior to computing said power spectrum. 
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6. The method of claim 2, wherein said analyzing each of said candidate 
frequencies to determine a best frequency by selecting n candidate frequencies 
associated with the n largest amplitude peaks in power spectrum which are not 
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harmonics of each other includes: 
bandpass filtering said segment of raw data; 

subsegmenting said bandpass filtered segment of raw data into individual heartbeat 
subsegments; 

computing an average heartbeat pulse from said individual heartbeat subsegments; 
evaluating said average heartbeat pulse against known plethysmographic pulse 

characteristics to obtain a quality measure for said average heartbeat pulse; 
repeating the above steps for each of said n candidate frequencies; and 
selecting a candidate frequency and average heartbeat pulse with highest quality 

measure to be said best frequency. 

7. The method of claim 6, wherein said bandpass filter is a finite impulse 
response (FIR) filter. 

8. The method of claim 7, wherein the coefficients of said FIR filter are 
adjusted to place the center frequency of said FIR filter at said candidate frequency. 

9. The method of claim 6, wherein said subsegmenting said bandpass 
filtered segment of raw data into individual heartbeat subsegments comprises: 
beginning each heartbeat subsegment at one quarter pulse width before a diastolic peak; 
ending each heartbeat subsegment at one quarter pulse width after the next diastolic 

peak, wherein each heartbeat subsegment overlaps with both a preceding 
heartbeat subsegment and a succeeding heartbeat subsegment, and each 
heartbeat pulse is centered within each heartbeat subsegment 

10. The method of claim 6, wherein said computing an average heartbeat 
pulse from said individual heartbeat subsegments comprises: 

calculating an average for each individual heartbeat subsegment; 
aligning corresponding individual heartbeat subsegment points from beginning to end 
for each of said individual heartbeat subsegments; and 
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calculating a modified average pulse by eliminating all data points outside x standard 
deviations from said average for each individual heartbeat subsegment to obtain 
said average heartbeat pulse. 

11. The method of claim 1 0, wherein x equals about two. 

12. The method of claim 6, further comprising minimizing dispersion of 
saturation values, wherein said minimizing dispersion of saturation values includes 
adjusting the phase relationship of red and IR pulses and calculating a minimum 
dispersion of saturation values. 

13. The method of claim 12, wherein adjusting the phase relationship of red 
and IR pulses includes calculating saturation values for a plurality of delays between 
said red and IR signals, wherein said plurality of delays are sequentially swept from 
about - 20 ms to about + 20 ms in 1 ms increments. 

14. The method of claim 12, wherein minimizing dispersion of saturation 
values includes: 

generating a histogram of saturation values with bin sizes of about y percent calculated . 
with said delays; 

locating a peak bin and its immediately surrounding two bins and calculating a measure 
of dispersion based on the percentage of saturation values which fall within said 
peak bin and it immediately surrounding z bins; and 

selecting an optimum delay corresponding to measure of dispersion with highest 
percentage. 

15. The method of claim 14, wherein y equals 0.25. 

16. The method of claim 14, wherein z equals two. 

17. A method of removing motion-induced noise artifacts from measured 
pulse oximetry signals prior to calculating functional pulsatile blood oxygen 
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concentration, Sp0 2 , comprising: 

acquiring a red data segment spanning a plurality of heartbeats from oldest to newest; 
frequency analyzing said red data segment for candidate red frequency components, one 
of which may be representative of a valid plethysmography pulse; 
5 analyzing each of said candidate red frequency components to determine a best red 
frequency component; 
outputting an average red pulse signal computed from said red data segment and said 

best red frequency component; 
repeating the above steps with said new red data segment; and 
10 substantially simultaneously performing the above steps with respect to an infrared data 
segment. 

18. A method of removing motion-induced noise artifacts from a single 
electrical signal representative of a pulse oximetry light signal comprising: 
15 (a) acquiring a segment of raw data spanning a plurality of heartbeats from said single 
electrical signal; 

(b) frequency analyzing said segment of raw data to determine two candidate 

frequencies with greatest amplitude, F, and F 2 , which are not harmonics of each 
other; 

20 (c) removing F 2 from said segment of raw data to obtain F 2 filtered data; 

~~(d) subsegmenting said F 2 filtered data into a plurality of subsegments each containing a 
single heartbeat pulse; 
(e) calculating a first average pulse signal based on averaging said plurality of F 2 
filtered subsegments; 
25 (f) removing F, from said segment of raw data to obtain F, filtered data; 

(g) subsegmenting said F, filtered data into a plurality of subsegments each containing a 

single heartbeat pulse; 

(h) calculating a second average pulse signal based on averaging said plurality of F, 

filtered subsegments; 

30 (i) determining which of said first average pulse signal and said second average pulse 
signal represents a valid plethysmography pulse; 
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(j) outputting said first average pulse signal or said second average pulse signal which 

represents a valid plethysmography pulse; and 
(k) repeating steps (b) through (j) with a new segment of raw data, 

5 19. The method of claim 1 8, wherein after step (b), frequencies F, and F 2 are 

compared with a valid frequency, Fo, from a previous signal processing cycle and F 0 is 
substituted for F 2 if neither Fj or F 2 are ± 0.1 Hertz of F 0 . 

20. The method of claim 18, wherein removing F 2 from said segment of raw 
10 data is accomplished by filtering said segment of raw data with a bandstop filter of a 

width about 0.5 Hertz centered at F 2 . 

21 . A method of removing motion-induced noise artifacts from a measured 
pulse oximetry signal comprising: 

1 5 acquiring a segment of pulse oximetry data from the measured pulse oximetry signal; 
filtering said segment of pulse oximetry data; 

performing frequency analysis on said filtered segment of pulse oximetry data to 

determine frequency components of said filtered segment of pulse oximetry 
data; 

20 selecting a frequency component with a largest power as an initial candidate heart rate 
response; _ 

subsegmenting said filtered segment of pulse oximetry data into a plurality of 

subsegments, wherein each subsegment represents an individual heartbeat; 

averaging data values for each subsegment at identical points along said length of each 
25 subsegment to form an average subsegment; 

determining whether said average subsegment represents a valid pulse oximetry signal; 

outputting valid average subsegments; 

updating said segment of pulse oximetry data; and 

repeating the above steps. 



30 



22. The method of claim 2 1 , further comprising measuring dispersion of 
saturation of latest pulse subsegment to determine whether said averaging step may be 
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skipped. 

23 . The method of claim 2 1 ; further comprising minimizing dispersion of 
saturation values prior to said outputting valid average subsegments. 

24. A circuit card for use in a pulse oximetry system to remove motion- 
induced noise artifacts from a measured pulse oximetry signal, said circuit card 
comprising: 

a circuit board for mounting electronic circuitry and interfacing with the pulse oximetry 
system; 

a processor mounted on said circuit board for processing at least one input signal 

according to instructions; and 
memory for storing a computer program, wherein said memory is operably coupled to 

said processor, and wherein said computer program includes instructions for 

implementing a method of removing motion artifacts from said measured pulse 

oximetry signal, said method comprising: 

receiving a segment of raw data spanning a plurality of heartbeats from said 

single electrical signal; 
analyzing said segment of raw data for candidate frequencies one of which said 

candidate frequencies may be representative of a valid plethysmography 

pulse; 

analyzing each of said candidate frequencies to determine a best frequency 

based on highest quality measure; 
outputting an average pulse signal computed from said segment of raw data and 

said best frequency; and 
repeating the above steps with a new segment of raw data. 

25. The circuit card of claim 24, wherein said processor is a digital signal 
processor. 



26. A pulse oximetry system for removing motion-induced noise artifacts 
from measured pulse oximetry signals comprising an input device, an output device, 



WO 00/22408 PCT/US99/23803 

-25- 

and a motion artifact circuit card, wherein said motion artifact circuit card comprises: 
a circuit board configured for coupling to said pulse oximetry system and 

communicating with other electronics within said pulse oximetry system; 
a processor mounted on said circuit board for processing at least one input signal 
5 according to instructions; and 

memory operably coupled to said processor for storing a computer program, wherein 
said computer program includes instructions for implementing a method of 
removing motion artifacts from said measured pulse oximetry signal, wherein 
said method comprises: 
1 o receiving a segment of raw data spanning a plurality of heartbeats from said 

single electrical signal; 
analyzing said segment of raw data for candidate frequencies one of which said 
candidate frequencies may be representative of a valid plethysmographic 
pulse; 

15 analyzing each of said candidate frequencies to determine a best frequency 

based on highest quality measure; 
outputting an average pulse signal computed from said segment of raw data and 

said best frequency; and 
repeating the above steps with a new segment of raw data. 

20 

27, The pulse oximetry system of claim 26, wherein said processor is a 
digital signal processor. 

28. A pulse oximetry system for removing motion-induced noise artifacts 
25 from measured pulse oximetry signals comprising an input device, an output device, 

and a motion artifact circuitry, wherein said motion artifact circuitry includes: 
a processor for processing at least one input signal according to instructions; and 
memory operably coupled to said processor for storing a computer program, wherein 

said computer program includes instructions for implementing a method of 
30 removing motion artifacts from said measured pulse oximetry signal, wherein 

said method comprises: 

receiving a segment of raw data spanning a plurality of heartbeats from said 
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single electrical signal; 
analyzing said segment of raw data for candidate frequencies one of which said 
candidate frequencies may be representative of a valid plethysmography 
pulse; 

analyzing each of said candidate frequencies to determine a best frequency 

based on highest quality measure; 
outputting an average pulse signal computed from said segment of raw data and 

said best frequency; and 
repeating the above steps with a new segment of raw data. 

29. The pulse oximetry system of claim 28, wherein said processor is a 
digital signal processor. 



WO 00/22408 



PCT/US99/23803 



1/11 



100 



110 



120 



130 



140 



Q Start ^ 



Acquire a Segment 
of Raw Data 



I 



Analyze Frequency 
Components 



I 



Determine Best 
Frequency 
Component 



I 



Output Average 
Pulse Signal 





New 
Segment \ 
fRawDatay^ 
? 



Yes 



End 



Fig. 1 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 



PCT/US99/23803 



2/11 



3.60E+05 




3.10E+05 J 

10.0 12.0 14.0 16.0 18.0 20.0 

Time (Seconds) 

Fig. 2 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 



PCT/US99/23803 



3/11 




BNSOOCJD: <W»_C0e2408A3J8> 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 



PCTAJS99/23803 



4/11 



3.6E+05 




3.1E+05 J 

11.0 13.0 1SJ0 17.0 19.0 

Time (Seconds) 
Raw IR Data Filtered 1R Data 

Fig. 4 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 PCTAJS99/23803 



5/11 



1.5E+04 




Data Points 

Fig. 5 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 



PCT/US99/23803 



6/11 



6000 




-6OO0 J 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Time (Seconds) 

Fig. 6 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 



PCT/US99/23803 



7/11 




700000 



■ i i j 1 i 1 i i *i 

30.00 31.00 32.00 33.00 34.00 35.00 36.00 37.00 38.00 39.0040.0040.80 

Time (Seconds) 



Fig. 7 



eNsooaa <wo_ooewoeA3je> 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 



PCT/US99/23803 



8/11 



10 



11 



12 



I/O 
Device 



Processor 



14 



Memory 



Fig. 8 



16 



18 



20 



Input 
Device 



Output 
Device 



i 



10 



I 



Motion Artifact Circuit Card 



Fig. 9 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 PCTAJS99/23803 



9/11 



22 



18 



20 



Input 
Device 



z 



Output 
Device 



%7 X 



Processor 



24 



Storage 
Device 



Fig. 10 



a«soocn>: <wo_ooe2«ieA3jB> 



SUBSTITUTE SHEET (RULE 26) 



PCT/US99/23803 



10/11 




Bandpass Filter 
Data Segment 
with FIR Filter 



I 



— © 



Subsegment 
Data 



Measure Dispersion 
of Saturation Based 
on Newest Pulse 




X 



720 



Is the 

Dispersion^. 
Low 
? 



Yes 



Compute Average 
Pulse 




Use Newest Pulse 

Instead of 
Computing Average 
Pulse 


<aU 







250 



Fig. 11 



SUBSTITUTE SHEET (RULE 26) 



WO 00/22408 



PCT/US99/23803 



11/11 



260 



270 



280 



290- 




Evaluate Pulse 
Shape Against 
Rules to Obtain 
Quality Measure 



I 



Minimize the 
Dispersion of 
Saturation Values 




Last 
Candidate 
Frequency 



Select Candidate 
Frequency with Best 
Quality Measure 



c 



Stop 



120 



Fig. 12 



SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



International application No. 
PCT/US99/23W3 



A. CLASSIFICATION OF SUBJECT MATTER 
IPC(7) :A61B 5/00 

US CL :600/336 

According to International Patent Classification (IPC) or to both national classification and IPC 

B. FIELDS SEARCHED 

Minimum documentation searched (classification system followed by classification symbols) 
U.S. : 600/310,322,323,330,336 



Documentation searched other than minimum documentation to the extent that such documents are included in the (itlds searched 



Electronic data base consulted during the international search '(name of data base and, where practicable, search terms used) 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



A.P 



US 5,490,505 A (DIAB et al.) 13 February 1996, the entire 
document. 

US 5,853,364 A (BAKER, JR. et al.) 29 December 1998, entire 
document. 



1-29 
1-29 



~] Further documents are listed in the continuation of Box C. | | See patent family annex. 



Special categories of cited documents: 

document defining the general state of the an which is not considered 
io be of panicular relevance 

earlier document published on or after the international filing date 

document which may throw doubts on priority claim's) or winch is 
cited to establish the publication date of another citation or other 
special reason fas specified) 

document referring to an oral disclosure, use. exhibition or other means 

document published prior to the international Cling date but later than 
the priority date csaaated 



T later document published after the international filing date or 

date and not in conflia with the application but cited to 
principle or theory underlying the invention 

"X" document of panicular relevance; the claimed invention rt—ot 

considered novel or cannot be considered to involve an inventive i 
when the document is taken a l on e 

"Y* document of panicular relevance; the claimed invention ci—tw I 

considered to involve ao inventive step when the 
combined with one or more other such documents, audi comb 
being obvious to a person skilled in rhe art 

~&~ document member of the same patent family 



Date of the actual completion of the international search 
28 FEBRUARY 2000 



Name and mailing address of the ISA/US 
Commissioner of Patents and Trademarks 
Box PCT 

Washington. D.C. 20231 
Facsimile No. (703) 305-3230 



Date of mailing of the international search report 



AuJfrfized 




ERIC F. WINAKUR 
Telephone No. (703) 308-3940 



